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Measurements of the azimuthal anisotropy of high-px neutral pion (ir°) production in Au+Au 
collisions at \fs NN = 200 GeV by the PHENIX experiment are presented. The data included in this 
paper were collected during the 2004 RHIC running period and represent approximately an order 
of magnitude increase in the number of analyzed events relative to previously published results. 
Azimuthal angle distributions of 7r°s detected in the PHENIX electromagnetic calorimeters are 
measured relative to the reaction plane determined event-by-event using the forward and backward 
beam-beam counters. Amplitudes of the second Fourier component (V2) of the angular distributions 
are presented as a function of tv° transverse momentum (pt) for different bins in collision centrality. 
Measured reaction plane dependent tv° yields are used to determine the azimuthal dependence of 
the 7T° suppression as a function of pt, 7?aa(A<^,pt). A jet-quenching motivated geometric analysis 
is presented that attempts to simultaneously describe the centrality dependence and reaction plane 
angle dependence of the tv° suppression in terms of the path lengths of hypothetical parent partons 
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in the medium. This set of results allows for a detailed examination of the influence of geometry in 
the collision region, and of the interplay between collective flow and jet-quenching effects along the 
azimuthal axis. 

PACS numbers: 21.65.Qr,25.75.-q,25.75.Dw 



I. INTRODUCTION 



Over the past few years, experiments at the Relativis- 
tic Heavy Ion Collider (RHIC) have established that a 
dense partonic medium_is formed in Au + Au collisions 

This medium thermalizes 
I S3) Eliis extremely opaque 



at y /I A 7^=200 GeV 
very quickly [l|, H, H,] 

to the passage of high-p T particles [I2|,|l3j], and the strong 
coupling of matter in the medium produces a system for 
which the ratio of shear visc osity to entropy (rj/s) ap- 
proaches zero 0, [H, [U, E3, El]- Much of the current 
focus is on the extraction of key transport and thermo- 
dynamic characteristics of the matter produced in these 
collisions. Measurements of high-px parton propagation 
in the medium as well as medium-induced modification 
of the fragmentation parton spectrum and its products 
provide a critical tool for probing medium properties. 

One of the most striking early results from RHIC 
was the observation of strongly suppressed production 
of high-px particles in central Au+Au events compared 
to appropriately scaled p + p collisions [HI, [l3| . High- 
px partons are formed from hard scattering between the 
initial colliding partons, and these partons fragment into 
two or more jets of hadrons. When propagating through 
a dense volume of deconfined matter, these high-px par- 
tons are expected to scatter from color charges in the 
medium, losing energy through a combination of gluon 
bremsstrahlung radiation and collisional energy transfer 
to partons in the medium. These radiated gluons even- 
tually fragment into hadrons at lower px, resulting in a 
depletion of the observed yields of hadrons at higher px- 

A useful way to quantify the suppression of high-px 
hadrons is the nuclear modification factor (-Z?aa) where 
the p + p cross section is scaled by the thickness function 
(Taa) of the two Au nuclei 



-Raa(pt) 



(l/NZ£)<PN?Jdprdy 
(T A a) x d 2 a^"/dp T dy ' 



PHENIX has measured a tt° Raa close to unity in both 
peripheral Au + Au collisions and Au collisions H(| , 
consistent with the expectation that these collisions 
would not produce an extended, dense medium. As the 
collisions become more central, Raa decreases to about 
0.2, indicating a stronger parton energy loss. Further- 
more, the measured ir° Raa is nearly constant as a func- 
tion of pt, for px ^ 5 GeV/c up to the highest currently 
accessibly p T , 20 GeVc 



'Deceased 
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These data can be well reproduced by models that cal- 
culate the energy lost by the hard scattered partons as 
they traverse the dense medium. The amount of energy- 
loss depends on the density of the medium [2l[ , so mea- 
surements of high-px hadron suppression provide con- 
straints on the transport coefficient (q 2 ), a measure of 
mean transverse momentum squared ( fc T ) transferred by 
the medium to a high-energy parton. However, multiple 
models with different physical assumptions can reproduce 
the measured Raa(pt) [12, The different models 

vary widely in how they include the crucial interference 
terms between multiple-scattering centers as well as the 
interplay between inelastic, elastic and flavor-changing 
processes during the parton's passage. 

To discriminate between these models we need to in- 
crease our experimental control of the path length, since 
the amount of energy lost by a high-px parton strongly 
increases with the distance traveled through the medium. 
A quadratic dependence on the path length is predicted 
for a static medium if the dominant energy-loss mecha- 
nism is the bremsstrahlung radiation of gluons surviving 
the destructive interference caused by multiple scatter- 
ing [i^, • For an expanding plasma the quadratic in- 
crease should be moderated to a linear dependence [24j ■ 

The centrality dependence of Raa (pt ) offers a probe 
of the path-length dependence of partonic energy loss. 
However, we can better test the path-length dependence 
by studying the azimuthal variation of the high-px sup- 
pression at a fixed centrality. Since the collision zone 
has a nearly elliptical shape in the transverse plane due 
to the non-central overlap of the colliding Au nuclei, par- 
tons that travel along the short axis of the nuclear overlap 
region lose less energy and should therefore be less sup- 
pressed. The key observable is then the two-dimensional 
modification factor i?AA(A(/>,px), where A</> is the an- 
gle of emission with respect to the event plane. The 
azimuthal dependence of the spectra can be also param- 
eterized by a Fourier expansion, where up to second order 
= Nq[1 + 2v 2 cos(2A<£)], with v 2 being called elliptic 
flow coefficient. While both quantities characterize az- 
imuthal asymmetries, historically and conceptually they 
have different roots. The notion of elliptic flow is pri- 
marily tied to lower px phenomena ("soft physics"), the 
domain where particle production is proportional to the 
number of participating nucleons (iV par t), and positive 
i>2 arises from the boost to the mean px in the direction 
where the pressure gradient is highest (along the reac- 
tion plane). Conversely, -Raa(pt) and i?AA(A0,px) are 
commonly used to describe high px behavior (hard scat- 
tering, which scales with the number of binary collisions 
-/Vcoll)- When Raa deviates from unity at high px, it be- 
comes a valuable probe of the loss of energy/momentum 
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in a particular direction. However, there is no clear sepa- 
ration between soft and hard regions, and both Raa an d 
V2 are well-defined in the entire momentum range, so in 
this sense V2 is sensitive to differential energy loss at high 
Pt- 

PHENIX has measured high-px v 2 for 7r° particles from 
Au+Au collisions [25[. The energy-loss models that re- 
produce Raa (Pt ) diverge in their predictions of the az- 
imuthal anisotropy at high px. They generally under- 
predict the observed azimuthal variation of .Z?aa(A0), or 
equivalently, are unable to describe the px dependence 
of i>2 over the full range of px where one would naively 
expect them to be applicable [26|, [27], Us}. These mod- 
els include the hydrodynamical evolution of the medium, 
and therefore the high-px probe loses energy in a medium 
that is becoming spatially isotropic with time. Several 
early papers noticed that the measured i>2 values were 
larger than what one would expect from a completely 
opaque almond-shape collision zone [1^, [3(| • Other early 
energy-loss calculations came close to reproducing the 
measured i?AA(A0) [H], H2], but in these the plasma ex- 
pansion was not taken into account, which resulted in 
unrealistically strong azimuthal anisotropy. Another cal- 
culation [33j] has reproduced Raa(^), but in this model 
the Au nuclei were parameterized as hard-spheres instead 
of using a more realistic Woods-Saxon density profile, 
and this mechanism artificially increases the azimuthal 
dependence of the energy density. 

One potential resolution of the problem with energy 
loss calculations not reproducing the measured azimuthal 
dependence of yields is a recent calculation that al- 
lowed the high-px parton to resonantly scatter with the 
medium [34( (and references therein), increasing the en- 
ergy lost by a parton at plasma densities that correspond 
to temperatures near the critical temperature. This pro- 
duces a sharper dependence of the energy-loss on the 
spatial variation of the medium's energy density and 
hence the model is able to simultaneously reproduce both 
-Raa(pt) and J?aa(A0). A critical check will be to exam- 
ine whether the same parameters work for the full range 
of collision centralities. 

In order to discriminate among all the models that at- 
tempt to reproduce i?AA(A0,px), the experimental chal- 
lenge is to extend the range and increase the precision of 
observations which can be used to test different energy- 
loss models. In this paper we extend the range of pub- 
lished data on 1?aa(A<^>) [25[ by a) reaching higher px, 
and thereby moving to a px region that is completely 
dominated by the fragmentation of hard partons and 
reducing the possible contribution of particles from re- 
combination [35j| . b) using finer bins in centrality, thus 
achieving less averaging of the path length, and c) reduc- 
ing the statistical and systematic uncertainties to further 
constrain models. We present in this article measure- 
ments using data collected during the 2004 RHIC run- 
ning period. These data represent a high-statistics sam- 
ple of Au + Au collisions (approximately 50 times that 
of the 2002 RHIC running period) and therefore extend 



our ability to measure -Raa(A0) and vi to much higher 
Px with better precision. 



II. EXPERIMENTAL DETAILS 

The data presented in this paper were taken by the 
PHENIX experiment [H in 2004 (RHIC Run-4), and 
represent the analysis of 821M minimum bias Au+Au 
collisions at ^/JWn — 200 GeV. The detectors involved 
in this analysis are the beam-beam counters 37] (BBC; 
triggering, centrality and reaction plane determination), 
the zero-degree calorimeter [38| (ZDC, centrality deter- 
mination) and the electromagnetic calorimeter [39] (EM- 
Cal, 7T° measurement). 

The BBCs are two groups of 64 hexagonal quartz 
Cerenkov radiator counters with photomultiplier readout 
surrounding the beampipe 144 cm up- and downstream 
("North" and "South") from the center of the nominal 
collision diamond, covering the 3 < \r)\ < 3.9 pseudora- 
pidity range and the full azimuth. Coincidence of signals 
in at least two photomultiplier tubes in both BBCs served 
as a minimum bias trigger and according to simulations 
it captured 92% of all inelastic collisions. The size of the 
total signal in the BBCs increases monotonically with col- 
lision centrality at this ^/snn- The collision vertex z was 
calculated from the difference between the fastest timing 
signals in the North and South BBCs, respectively, with 
a < 2.0 cm resolution. Only events with \z\ < 30 cm 
were analyzed. 

The ZDCs are small tungsten/scintillator hadron 
calorimeters with quartz fiber lightguides and photomul- 
tiplier readout, located between the beampipes at 18 m 
North and South from the collision point. They measure 
non-interacting "spectator" neutrons in a cone of about 
2 mrad, and their signal is double-valued as a function 
of centrality (it is low in very central and very peripheral 
collisions but large at mid-centrality) . The correlation of 
ZDC vs BBC signals resolves this ambiguity and allows 
a precise measurement of the true centrality for all but 
the most peripheral collisions. 

The reaction plane (spanned by the beam direction and 
the impact vector of the colliding nuclei) is determined 
event-by-event from the azimuthal charge distribution in 
the BBCs, after taking into account small nonuniformi- 
ties (in the response of individual radiators, PMTs, elec- 
tronics, etc.), using the assumption that over a large 
number of events the <fi distribution of per-event reac- 
tion planes should be uniform. Due to the large rapidity 
gap between the central arm (I77I < 0.35) where 7r°s are 
measured, and the BBCs where the reaction plane is es- 
tablished, we assume that the reaction plane is unbiased 
and free from auto-correlations. However, the relatively 
coarse granularity of BBCs affects the resolution. Note 
that in this analysis precise knowledge of the reaction 
plane resolution, which depends strongly on centrality, 
is crucial. This will be discussed in detail in the next 
Section. 
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Neutral pions are measured by reconstructing their de- 
cay photons (tt° -> 77) in the EMCal. The EMCal 
consists of 8 sectors at midrapidity (\r]\ < 0.35), cov- 
ering a total of 2 x 90° in azimuth. Six sectors are 
lead/scintillator (PbSc) sampling calorimeter with pho- 
tomultiplier readout and 5.5 x 5.5 cm 2 granularity, two 
sectors are lead/glass (PbGl) Cerenkov counters with 
4x4 cm 2 granularity and photomultiplier readout. The 
two detectors are 18X and 16X radiation lengths deep, 
respectively, both ensuring essentially full containment 
of electromagnetic showers in the relevant energy range. 
The in situ energy resolution is well reproduced by sim- 
ulation both in PbSc and PbGl: the ir° peak positions 
and the widths both agree with the data to better than 
1 MeV over the entire momentum range. Therefore, the 
error on the energy (and momentum) scale is less than 
1%. Timing resolution er t is ~ 450 ps and ~ 650 ps 
for the PbSc and PbGl, respectively, allowing the rejec- 
tion of neutrons and antineutrons up to a few GeV/c 
transverse momentum, which would otherwise be a ma- 
jor source of neutral showers up to a few GeV energy. 
At sufficiently high transverse momenta, decay photons 
from a nearly symmetric (E 7l sa E l2 ) decay may produce 
showers in the calorimeter that start to merge into one 
reconstructed cluster. In the PbSc this effect is first vis- 
ible around pt ~ 10 GeV/c, at the upper end of the pt 
region considered in this paper. Due to its higher gran- 
ularity and smaller Moliere-radius the PbGl is immune 
to this "merging" problem up to px ~15 GeV/c. The 
hadronic response, timing properties and other sources of 
systematic errors are very different for the two calorime- 
ter types. Therefore, when extracting the ^-integrated 
i?AA, which serves as absolute normalization, the PbSc 
and PbGl were analyzed separately and the results com- 
bined to decrease the total systematic uncertainty. 



III. DATA ANALYSIS 

A. Centrality 

As mentioned, the minimum-bias trigger in the Run-4 
PHENIX configuration is supplied by the BBCs, and the 
correlation of the charge deposited in the BBCs with en- 
ergy deposited in the ZDCs provided a determination of 
the centrality of the collision. The elliptic flow measure- 
ment presented in this paper is measured in seven bins 
of the centrality range 0-92%, with lowest correspond- 
ing to the most central: 0-5%, 5-10%, 10-20%, 20-30%, 
30-40%, 40-50%, and 50-60%. In addition, the combined 
ranges 0-20%, 20-40%, 40-60%, and minimum bias bins 
are included. For the yields with respect to the reaction 
plane, the centralities presented are 0-10%, 10-20%, 20- 
30%, 30-40%, 40-50%, 50-60%. Finally, the R AA versus 
nuclear path length result excludes the most central bin 
due to its smaller intrinsic ellipticity (the average path 
length is insensitive to A0). 
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FIG. 1: Reaction plane resolution correction as a function of 
centrality. 



B. Reaction plane determination 

The technique used to determine the reaction plane 
on an event-by-event basis is the same used in previous 
PHENIX analyses 0, 0, E|- The quartz radiators of 
each counter are arranged in approximately concentric 
circles around the beam axis. The light collected in the 
photo-multiplier tubes (PMTs) allows for an estimate of 
the number of charged particles passing through the de- 
tector. 

The number of charged particles at a given PMT po- 
sition, Ni, is weighted in a manner to reduce the bias of 
the inner rings and used to measure the orientation of 
the reaction plane from the formula 



tan(2*) 



J2i WjNj sin(20j) - (J2j WjNj sin(2<^)) 
£\ w z N, cos(24) - (£\ Wi Ni cos(2<^)) ' 



(1) 

where <pi is the nominal azimuth of the radiator. Subtrac- 
tion of the average centroid removes biases due to various 
detector effects. A final flattening technique is used to 
remove the residual non-uniformities in the distribution 
of angles. 

To estimate the resolution of the reaction plane mea- 
surement, we use the sub-event technique [42|]. The ap- 
proach consists of dividing the event up into two sub- 
events roughly equal in size. The two individual BBCs 
provide a natural sub-event division, so we analyze the 
distribution of event-by-event differences between the re- 
action plane angles measured in the north and south 
counters, Afy = — ^s- In t he p resence of pure flow, 
this distribution takes the form 14211: 
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dN 
dA^ 



-(1 + x 2 ) + z [I (z) + Lo(z)] + x 2 + L x (z)] 



(2) 



where z = x 2 cos (2A^) and the functions J„ and L n are the modified Bessel functions of the first kind and modified 
Struve functions, respectively. The parameter x describes the dispersion of the flow Q vector and thus determines 
the correction required for the reaction plane resolution. Since A "J represents the whole-event difference distribution 
and we are dealing with sub-events with roughly half the multiplicity of the event, we replace x ^ x/V% m Eq. [2] 
and fit this function to the measured A\& distribution to extract y. The resulting value is then used to evaluate the 
resolution of the event-plane of n th order (42j: 



In-l 



(3) 



where the true reaction plane orientation is denoted 
by *&rp and the observed orientation by "J. Figure Q] 
shows the resolution correction obtained using the above- 
described procedure as a function of centrality. Both 5% 
and 10% wide bins are shown for comparison. 

Eq. [2] is derived under the assumption that the az- 
imuthal distributions are free of non-flow effects. Due to 
the large rapidity gap between the BBCs and the cen- 
tral arm region, it is expected that particles observed in 
the BBCs have no correlation with those measured in 
the central arm detectors. PYTHIA [43| studies have been 
used to confirm that jets observed in the central arm have 
negligible effect on the reaction plane measurement from 
the BBCs 03. 



C. Neutral pion measurement 

Measurement of neutral pions has played a critical role 
in the study of highjgr phenomena at RHIC, and espe- 
cially by PHENIX pj, El Sj| . The two-particle decay 
channel ir° — > 7 + 7 provides a clean signal of identified 
hadrons out to the highest pt regions. 

EMCal showers are found by clustering contiguous 
towers with energy above a threshold energy (10 MeV) 
and requiring at least 50 MeV in the tower with highest 
energy deposit. The impact position is calculated from 
the positions of the participating towers weighted by the 
logarithm of deposited energy. The energy of the cluster 
is corrected for non-perpendicular incidence - the angle 
being derived by assuming a straight path between the 
actual vertex and the calculated impact point - as well 
as nonlinearities [25j . In high- multiplicity events such as 
central Au + Au collisions, there is an increasing proba- 
bility for clusters to overlap (one tower accumulates en- 
ergy from more than one particle), which can distort an 
energy measurement from a simple sum over contiguous 
towers. To mitigate this effect, the EMCal clustering al- 
gorithm also provides a quantity called ecore, which is de- 
termined by extrapolating the "core" energy represented 
by the central four or five towers in the cluster, assum- 



ing an electromagnetic shower profile. The energy- and 
impact angle- dependent shower profile is a model devel- 
oped from and checked against beam test data. In this 
way, ecore provides a more realistic measurement of the 
shower energy, less prone to contributions from acciden- 
tal overlaps (particles hitting close enough they deposit 
energy in the same towers) than a simple energy sum of 
participating towers would be. We use this ecore for the 
energy of reconstructed clusters in this analysis. 

The invariant mass of a photon pair 7$ , jj as measured 
in the EMCal is calculated from the energy of the clusters 
and their measured position: 

™7,7, = ^(K+ P V = ^/^COMl ».,) (4) 

where 9ij is the opening angle between the two photons 
and m 7i7 - is equal to the ir mass for photons from the 
decay of the same ir°. Since the photons from the ir° are 
not tagged, such pairs have to be formed from each pho- 
ton pair in the event where the pair momentum falls in a 
particular px bin, and some of these pairs might acciden- 
tally reproduce the ir° mass as well (combinatorial back- 
ground), particularly at lower p^ and higher centralities 
(multiplicities). Since 7r°s cannot be uniquely identified, 
raw 7T° yields are extracted statistically, by subtracting 
the combinatorial background from the invariant mass 
distribution. 

A well-known technique to reduce the combinatorial 
background is to place a cut on the energy asymmetry of 
the pair, as defined by: 

a =^fi=/3|cosn (5) 

Because the angular distribution da /d cos 8* of the pairs 
in the rest frame of the 7r° is uniform, the asymme- 
try distribution should be flat. However, due to the 
steeply falling photon spectrum, fake (non-correlated) 
pairs which still give the proper tt° mass are strongly 
peaked towards a = 1. A pair of clusters in the EMCal 
is considered a neutral pion candidate only if the pair's 
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asymmetry is less than 0.8. In addition, the two pho- 
tons are required to be separated by at least 8 cm for the 
combination to be considered as a 7r° candidate. 

There remains a non-trivial background contribution 
which passes these cuts: pairs of photons from different 
7r°s, or, more generally, from pairs of uncorrelated clus- 
ters which pass the photon identification cuts and acci- 
dentally give an invariant mass near the true 7r° mass. 
This remaining combinatorial background is estimated 
and subtracted using the event mixing method. The 
procedure involves forming pairs from different events, 
which will by definition be uncorrelated. Each photon 
candidate is combined with all the photon candidates in 
previous events stored in memory. In order to replicate 
the background from uncorrelated pairs within the same 
event as closely as possible in the mixed events, mixing is 
performed within bins of ccntrality, vertex z position, and 
reaction plane orientation. Since all events analyzed are 
minimum-bias, no special steps are needed to avoid the 
distortions of the mixed-event background by the trig- 
ger requirement. All cuts applied to the combinations of 
same event pairs are also applied to mixed-event pairs. 
The number of events buffered determines the statistics 
of the event-mixed distributions, chosen as a tradeoff be- 
tween desired statistical accuracy and computational re- 
sources. The data presented in this article are mixed 
with five previous events (in each centrality, vertex, and 
reaction plane bin). 

For a given bin, the mixed-event mass distribu- 
tion is normalized to the same-event distribution in a 
region away from the ir° mass peak. The normalization 
region is 0.25-0.45 GeV/c 2 forp T < 6.0 GeV/c and 0.21- 
0.45 GeV/c 2 otherwise. Fig. [5] shows an example of this 
subtraction process for two pt ranges in two centrality 
bins. 

The scaled background distribution is then subtracted 
from the same-event pair distribution. The subtracted 
result thus represents a sample of real ir°s. The peak is fit 
to a Gaussian to determine its width and mean position. 
The raw yield of 7r°s is determined by integrating the 
counts in a window of ±2cr around the mean. The width 
and mean are recorded and parameterized as a function 
of pt and centrality based upon this ^-integrated, large 
sample. The positions and widths from this parameteri- 
zation arc then used when we extract the (much smaller) 
raw yields in bins of angle A(f> with respect to the reaction 
plane. The maximum variation of the yields (multiplici- 
ties) with Atfi is only about a factor of 2, and therefore the 
means and widths are not expected to change substan- 
tially. Furthermore, the statistics are much poorer in the 
A(f) bins, which would make individual ir° fits unreliable. 

There is a residual background in the invariant mass 
distributions even after the mixed-event distribution has 
been removed, especially at lower pt (below ~ 2 GeV/c). 
This is due to correlations that event mixing cannot re- 
produce, like the "sub-event structure" due to the pres- 
ence of jets or multiple, close-by showers from an annihi- 
lating anti-neutron, or imperfections of the reconstruc- 



tion algorithm, such as cluster merging, cluster split- 
ting, and a host of other contributions. Much of the 
residual background is excluded by starting the fit at 
0.09 GeV/c. What is left is accounted for by includ- 
ing a first-order polynomial in the fits to the (already 
background-subtracted) invariant mass distribution, and 
subtracting its integral from the raw 7r° yield (see Fig. [2]). 
In the more central events, the peak deviates slightly 
from gaussian on the high mass side, due to overlapping 
clusters. The use of ecore mitigates this effect, and the 
systematic uncertainty on yield extraction arising from 
the remaining asymmetry has been estimated to be 3- 

4% M- 

D. Elliptic flow measurement 

To obtain the azimuthal angle dependence of ir° pro- 
duction, we measure raw 7r° yields in a given bin 
as a function of the tt° angle with respect to the re- 
action plane orientation in six equally-spaced bins of 
A<fi = 0(vr°) — ^rp covering the range < A</> < ir/2. 
The 7T° yields are measured in each Acj) bin using the 
same procedure described in llll CI for the reaction-plane 
inclusive measurement except that the mass fits are not 
performed in each A(j) bin. Instead, the peak integration 
window is set ±2cr around the mean where the width and 
mean are taken from the inclusive analysis. The result- 
ing raw 7T° angular distribution dN/dA(f> can then be fit 
to determine the strength of the modulation in the yield. 
Because the PHENIX BBCs have uniform azimuthal cov- 
erage, the 7T° measurements have uniform acceptance in 
A(f> when averaged over a large event sample, despite 
the limited azimuthal acceptance of the PHENIX elec- 
tromagnetic calorimeters. 

Assuming elliptic flow is the dominant source of Acf> 
variation in the tt° yields, we perform a fit to the angular 
distributions of the form 

dN 

— =N (l + 2 V2 mcas cos2A^). (6) 

We use an analytic linear \ 2 fitting procedure that 
matches the integral of Eq. [6] over each of the A0 bins to 
the measured ir° yield within the corresponding bin. In 
the definition of the y 2 function we account for non-zero 
covariances between the yields in the different A<fi bins 
resulting from the limited acceptance of the calorime- 
ters. These covariances have been evaluated separately 
for each and centrality bin. Examples of the raw 
dN/dA(j) distributions and the results of the % 2 fits are 
shown in Fig. [3J The resulting t>2 lcas values are then cor- 
rected upward to account for reaction plane resolution 
using correction factors described in Section ITlI B I 
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FIG. 2: Invariant mass distributions at moderate pt and two different centralities. Left panels: same event and normalized 
mixed event distributions. Right panel: the subtracted distributions, which are then fitted with the sum of a first degree 
polynomial and a Gaussian. 
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E. Raa measurement 

The nuclear modification factor Raa has played a crit- 
ical role in understanding energy loss mechanisms. Raa 
is defined as 



Raa(j>i 



(l/N^)d 2 NZ A /dp T dy 
(T AA ) x d 2 a£/dp T dy 



(7) 



to the data represents an approximation to a full un- 
folding procedure that becomes exact when the distribu- 
tion is purely cos (2A<fi) in form. We have checked for 
a few cases that a full unfolding procedure applied to 
the measured dN/A(j) distributions using singular value 
decomposition regulation of the response matrix repro- 
duces the correction in Eq.[l2l From the corrected ratios, 
i? corr (A<^,p T ), we use Eq.CTOlto obtain #aa(A^,Pt)- 



where (Taa) is the mean Glauber overlap function for 
the centrality being analyzed: 



(T. 



A A 



/T AA (b)db 



/(I 



ci T AA (b) 



)db 



(8) 



from which the mean number of binary nucleon-nucleon 
collisions can be calculated, (N coll ) = <™ e ' (T, 



aa 



For each px bin, we can calculate the ratio 



i?(A</> 4 ,p T ) 



jV(A0 t ,p T ) 



(9) 



where N(A4>i,pr) is the number of 7r°s observed in the 
given (A0j,pr) bin. Since the BBC is azimuthally sym- 
metric the PHENIX acceptance has no Acf> dependence, 
there should be no azimuthal dependence to efficiency 
and acceptance corrections. As a result, 



i?AA(Ac^,pr) = R(A(j>i,p T ) x Raa (pt ) ■ 



(10) 



Thus, we can use measured inclusive Raa {pt ) to con- 
vert R(A(f>i,p-r) to i?AA(A0i,px)- Since the detector 
efficiency and acceptance corrections are already con- 
tained in -Raa(pt), there is no need to apply them to 

Prior to calculating i?AA(A</>i,px) we correct the ra- 
tios i?(A0i,px) for the finite reaction plane resolution 
using an approximate unfolding technique. For a pure 
flow A<j) distribution, we can express the influence of the 
resolution broadening on the measured A<p distribution 



i? meas (A<^,p T ) = i? tluc (A0. i ,p T ) 



1 



2u moas 



COs(2A0) 



2^ orr 



cos(2A0) 

(11) 

where according to the results of Section IIIIBI u 2 ncas = 
W2 Orr /( C0s2 ( v I / - ®rp))- Then, if the measured A<p dis- 
tribution resulted from pure elliptic flow, it could be cor- 
rected back to the true distribution by 



fl corr (A&,p T ) = i? roea8 (A^,p T ) 



1 



2u§ orr cos(2A0) 



1 



2w 2 mcas ■ 



! cos(2A0) 
(12) 

As shown above, the general features of the measured ir° 
A(j> distributions are well-described by pure cos (2A</>) 
modulation. However, we wish to preserve in our mea- 
surements of the azimuthal dependence of the tt° pro- 
duction the full shape of the measured A(f> distribution, 
including possible small non-elliptic contributions. For 
this purpose, the correction described in Eq. [T^] applied 



IV. RESULTS 
A. Elliptic flow coefficient 

The results of the V2 measurements using the methods 
described in Sec. MIDI are presented in Fig. [4] as a func- 
tion of j>t for different centrality bins. The data points in 
the figure are plotted at the mean ir px in bins of width 
Ap T = 0.5 GeV/c for p T < 4 GeV/c and Ap T = 1 GeV/c 
for pt > 4 GeV/c. The error bars shown on the vi data 
points were obtained by multiplying the raw vi fit errors 
(see Section MI Dl) by the same reaction plane resolution 
correction factor applied to the V2 values themselves. The 
error bars, then, represent uncorrelated statistical errors 
on the measured vi values arising from statistical errors 
on the dN/dA(f> data points used in the fits (these er- 
rors would be categorized as Type A uncertainties in the 
framework described in [45| or px-uncorrelated) . System- 
atic errors on the vi measurements due to the reaction 
plane determination procedure and from systematic un- 
certainties in the reaction plane resolution correction are 
represented in Figs. 14151 by filled boxes, which for most 
data points are similar in size or smaller than the data 
points (these uncertainties are classified as Type B [45[ 
or px-correlated) . 

Figure [5] shows U2(pt) for four centrality ranges, ob- 
tained by combining data from the centrality bins shown 
in Fig. [U The corrected dN / dA<f> distributions from indi- 
vidual centrality bins were summed over a given central- 
ity range and then fit to obtain the corrected %)% values 
shown in Fig. [5] The reaction plane resolution correc- 
tion produces correlated errors in the corrected dN/dAcj) 
distributions for each original centrality bin, and these 
correlated errors persist in the summed dN/dAcj) distri- 
bution. These correlated errors are not included in the 
statistical errors for the fit to the combined dN/ dAcp dis- 
tribution, but their impact on the final V2 value is es- 
timated separately by evaluating the changes in the i>2 
fit parameter that result from adding to the summed 
dN/dA<p values ±1 a of the correlated errors. Since this 
estimated uncertainty results from the statistical uncer- 
tainties on the v% values for the original centrality bins, 
we include the la bounds obtained from this procedure 
in the statistical error on the V2 values for the combined 
centrality bins. Systematic errors for the combined bins 
are plotted similarly to those in Fig. [4] 

The results presented here nearly double the px range 
of previous PHENIX ir° V2 measurements from RHIC 
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Run- 2 [40(| . Those measurements are shown for com- 
parison purposes in Fig. [5] Good agreement is seen be- 
tween the Run-4 measurements presented here and the 
Run-2 results except in the 40-60% centrality bin where 
the new vi measurements are systematically higher by 
~ 30%. This difference is attributed to improved reac- 
tion plane resolution corrections for the 40-60% central- 
ity bin resulting from the combining of corrected dN/Acf> 
distributions from smaller centrality bins. This summing 
procedure better handles the rapid variation of reaction 
plane resolution with centrality in mid-central to periph- 
eral collisons. Furthermore, we have cross-checked the 
procedure using 5% bins, verifying the combined result 
reproduces the data analyzed in wider bins. The previ- 
ous Run 2 analysis did not have a sufficiently large data 
sample to allow the use of separate 40-50% and 50-60% 
bins, and therefore the reaction plane resolution correc- 
tion was necessarily less accurate. The measured Vi val- 
ues presented in Fig.[5]and Fig. [5] are also consistent with 
previously published PHENIX charged pion V2 measure- 
ments pj,l40l|. 

The results in Figures [4] and [5] show a rapid increase of 
V2 with increasing px at low px, a maximum in the range 
2 < pt < 3 GeV/c, and then at higher px a decrease of i>2 
with increa sing »r. An increase in V2 at low pt is well- 
established [iy, [4l], Hl| and is understood to result from 
the collective elliptic flow of the medium generated by 
the initial spatial anisotropy of the collision zone. Hydro- 
dynamical models have been successful in quantitatively 
describing the pion V2(px) in the region px < 1-5 GeV/c. 
However, it has also been well-established that the pion 
V2{pt) deviates from the hydrodynamic prediction above 
1.5 GeV/c, a result that is understood to imply the con- 
tribution of hard processes, distortions of the spectrum 
due to recombination at freeze-out, and/or effects from 
dispersive hadronic evolution after freeze-out. Thus, a 
change in the variation of V2 with px near px ~ 2 GeV/c 
is not unexpected. If the large V2 values at lower px 
are interpreted as resulting from soft, collective mecha- 
nismsm, then a decrease in V2 for px > 3 GeV/c sug- 
gested by the data in Fig. [4] would naturally reflect an 
increasing contribution of hard processes with smaller V2- 

To statistically test the significance of the decrease of 
V2 with px, we show in Fig. [5] the results of linear fits to 
the high-px V2 values for the 20-30% and 30-40% central- 
hies. The panels on the left-hand side of Fig. [S] display 
a series of fits beginning at higher px values, the first fit 
starting at the px near the maximum v%, px = 2.5 GeV/c. 
The right-hand panels show the 1-a limits of the func- 
tions for the three fits in the series, calculated from the 
1-a variation of the two fit parameters (and including the 
covariance between them) . The results of the fits indicate 
that the decrease of V2 with px at higher px is statistically 
significant, though the data points for px > 5 GeV/c are 
not sufficient by themselves to establish a trend. We 
can state, however, that the points for px > 5 GeV/c 
are consistent with the linear decrease obtained includ- 
ing the lower px points. A question we would like to 



answer, then, is whether the data show any indications 
of devation from a monotonic decrease in i>2(px) indicat- 
ing the transition to a quenching-dominated azimuthal 
variation. 

A complete understanding of ^(px) over the measured 
px range therefore requires the treatment of the transi- 
tion from soft to hard dominated physics. According to 
the above discussion, in the px range where V2 is maxi- 
mum, particle production is likely not dominated by hard 
processes and the reduction of V2 with increasing px in- 
dicates increasing hard-scattering contributions (or de- 
creasing soft contamination). Motivated by this general 
argument, we have attempted to describe the results in 
Figs. 2] and [3] by a functional form 



v2i.P1) = 



(px/Af 



1 + (p T /A) r 



1 



3T • (13) 



The first term is intended to describe a rapidly rising and 
saturating soft V2 resulting from collective motion while 
the second term represents a rapidly falling soft/hard 
ratio. The additive constant in the second term repre- 
sents an asymptotic V2 that could describe a constant or 
slowly varying azimuthal-dependent quenching. We show 
in Figs. 0[8]the optimum fits to the full set of i>2(px) val- 
ues in the different centrality bins and the result of la 
variation of the fit parameters taking into account the 
complete covariance matrix from the fits. 

The fits to the data show that the measured px de- 
pendence of the 7T° V2 is qualitatively compatible with 
a description of the of low and intermediate px region 
in terms of a collective flow modulation diluted by a 
decreasing relative soft contribution with increasing px. 
Assuming this picture, it is then important to determine 
at what px the contamination from the soft production 
no longer dominates the measured A<fi variation of 7r° 
yield. For most of the centrality bins, the fits suggest 
that V2 decreases over most of the measured px range 
albeit with a decreasing slope at higher px. The cen- 
tral bins are compatible with V2 continuing to decrease 
beyond the measured px range although the la uncer- 
tainty bands also accommodate V2 saturating within the 
measured range. The more peripheral bins (30-40% and 
40-50%) suggest that the V2 has reached a nearly px in- 
dependent value by ~ 5 GeV/c. The 50-60% centrality 
bin has sufficient fluctuations that little can be inferred 
from the px dependence of V2 in that centrality bin. In 
all of the centrality bins, the data are consistent with a 
smooth reduction of t>2(px) from a maximum to a non- 
zero value at high px with that value increasing in more 
peripheral collisions as would be expected from quench- 
ing in an increasingly anisotropic medium. While the 
functional form in Eq. [T3] can describe the px variation 
of V2 within the range of the current data and within the 
statistical fluctuations of the data points, it is possible 
that this description will fail at higher px with improved 
statistics. In fact, a statistically significant deviation of 
t>2(px) from the form in Eq. 1131 might provide the most 
direct evidence of the dominance of quenching effects in 
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FIG. 6: (Color online) vi versus px for 20-30% and 30-40% centralities, with fits of the high-px data to a first order polynomial. 
From right to left, the first panel shows the series of fits, with each fit starting with a successively higher px. The second, third, 
and fourth panels show selected fits with uncertainty bands based on the 1-a variation of the fit parameters, including their 
covariance. 



determining V2- 



Nuclear modification factor with respect to the 
reaction plane 



The nuclear modification factor as a function of A</> 
for six centrality bins is shown in Figs. lTJl[l"4l The closed 
circles represent the A0-dependent measurements de- 
scribed in this paper while the open circles positioned at 
A<ft = 7r/4 represent the inclusive Raa measurement [l9| . 
In both cases statistical uncertainties (i.e. Type A) are 
represented by the error bars. For the inclusive Raa 
measurement, the total systematic uncertainties (or Type 
C [45]) are shown by the boxes. The upper and lower 
la ranges of the correlated statistical uncertainties (i.e. 
Type B) on the -Raa(A</>) measurements resulting from 
the reaction plane resolution correction are shown by the 
(blue) solid and (red) dashed lines. For all bins except the 
0-10% centrality bin a dotted line is plotted at Raa = 1 
for reference. We note that by construction, the average 
-Raa(A</>) from the reaction plane dependent measure- 
ments must be equal to the inclusive Raa- 

The results in the Figs. IMT41 show that the in-plane w° 
suppression is generally weaker and varies more rapidly 
with pt than the suppression for 7r°s produced at larger 
angles. As the collisions become more peripheral (for 
example, 50-60%), the small suppression seen in the in- 
clusive Raa almost vanishes for 7r°s emitted close to the 
reaction plane. In a previous analysis 25], it was ob- 
served that the in-plane Raa even exceeded unity for 
peripheral collisions; these data exhibit no such enhance- 
ment. However, the results presented in this article agree 
within systematic errors with previously reported data. 

The i?AA(A0) results are combined in Fig. [15] that 
shows the pr dependence of the Raa in each of the six 
bins included in this analysis. We use a semi-log scale 
so that the reduction of the A</>-integrated -Raa in more 
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FIG. 7: 7T° v 2 versus p T for centralities 0-5%, 5-10%, 0-10%, 10-20%, 20-30%, 30-40%, 40-50%, and 50-60%. The arrow in the 
50-60% panel shows the lower limit of the uncertainty on the data point, which lies outside the bounds of the plot. The solid 
lines represent the fit to the data, Eq. 1131 The dashed lines represent the 1 a deviations of the fit function. See text for more 
details. 



central collisions does not confuse the interpretation of 
the results. For clarity, the results from the 20-30%, 30- 
40%, 40-50%, and 50-60% centrality bins are shown on 
linear plots in Fig. [TBI 

The -Raa(pt) results exhibit a peak near 2 GeV/c, 
which becomes more prominent for more central colli- 
sions. The peak is strongest in the 0-10% bin where there 
is little modulation of the A<fi distributions at low or high 
Pt, so the peak cannot be directly attributed to elliptic 
flow. The peak in Raa near 2 GeV/c is much weaker 
in the more peripheral (40-50% and 50-60%) centrality 
bins, particularly for 7r°s produced at larger A(/>, and the 
primary variation seen in these peripheral bins with in- 
creasing A(j> is a reduction in Raa that is only weakly px 
dependent. 

For the intermediate centrality bins (10-20% through 
30-40%) the peaking in Raa is seen in all A<fi bins, but 



is much stronger in the in-plane bins. For these inter- 
mediate centralities and for values above the peak 
in Raa (pt ^ 3 GeV/c), the Raa for ir°s produced at 
angles normal to the reaction plane is nearly constant 
with pt while the Raa for 7r°s produced at small angles 
from the reaction plane decreases rapidly with increas- 
ing pt. The near constancy of the out-of-plane Raa to- 
gether with the rapid reduction in in-plane Raa indicates 
that in the intermediate centrality bins, the 1)% and in- 
clusive Raa decrease simultaneously with increasing 
such that R aa(tt /2,px) = -Raa(pt)(1 — 2W2) is approxi- 
mately constant. We will argue below that a correlation 
between Raa and i>2 may naturally result from the un- 
derlying physics responsible for the azimuthal variation 
of the particle yields. However, we observe that a simul- 
taneous reduction in integrated Raa and V2 suggested 
by the more central Raa (pt ) data would be contrary to 
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naive energy loss expectations since smaller Raa would 
imply stronger quenching in the medium which would, in 
turn, imply larger variation between in-plane and out-of- 
plane quenching . 

A similar implicit correlation between integrated Raa 
and V2 is seen in the centrality dependence of the 
.Raa(A</>) results. These are re-plotted in Fig. [T7] as a 
function of iVpart for three A<fi bins - the bins closest 
to and further from the reaction plane and one of the 
intermediate bins. For N palt > 100, the out-of-plane 
Raa values are nearly independent of centrality while 
the in-plane Raa values decrease rapidly with increas- 
ing centrality. This result would have a natural geomet- 
ric interpretation for ir° production dominated by hard 
scattering and jet quenching. The length of the medium 
normal to the reaction plane varies only slowly with cen- 
trality except in the most peripheral collisions. Then, if 
the 7T° suppression is determined primarily by the path 
length of its parent parton in the medium, the ir° Raa 
would be nearly constant. Following the same argument, 
the yield for pions in the direction of the reaction plane 
would be much less suppressed in non-central collisions 
due to the short path lengths of the parent partons in 
the medium. However, with increasing centrality and de- 
creasing anisotropy of the collision zone, the in-plane par- 
ton path lengths would grow to match those in the out- 



of-plane direction. Thus, if the ir° suppression depended 
primarily on path length, the in-plane Raa would natu- 
rally drop to match the out-of-plane values reproducing 
the behavior of Fig. [17] In order to better see the dif- 
ference between the in- and out-of-plane behaviors, these 
data are also plotted on Fig. [T5] with a semi-log scale. 

One difficulty with this geometric interpretation of the 
.Raa(A0) results given above is that the trend in the data 
that it is supposed to explain persists down to low px, 
where the V2 values are too large to be accounted for via 
perturbative or formation time based energy loss scenar- 
ios [n m in m, m m. That ^ coupled with the 

pronounced peaking in Raa{^<P) near 2 GeV/c suggests 
that physics other than hard scattering and jet quenching 
must be invoked to explain the ir° yields at intermediate 
Pt. However, the fact that the out-of-plane yields show 
less pronounced peaking near 2 GeV/c, that they vary 
little as a function of pt above 3 GeV/c, and that they 
vary little with centrality for N part > 100 could be inter- 
preted to imply that the ir° suppression at angles normal 
to the reaction plane more directly represents the effects 
of quenching of hard quarks and gluons while the yield 
of 7r°s produced more closely aligned with the reaction 
plane is enhanced by the collective motion of the system. 
That additional enhancement could either be due to soft 
hadrons being boosted to larger px values by the collec- 
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FIG. 12: 7T° Raa versus angle of emission with respect to the reaction plane for 30-40% centrality. Colors/data points as in 
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tive elliptic flow or could result from weaker quenchin g fo r 
partons moving in the direction of the flow field [13, • 
Simultaneous description of the in-plane and out-of-plane 
behavior is a sensitive test of energy loss models. 

The v 2 values presented in Figs. 05] also peak near 
2 GeV/ c, but the locations of the maxima in v 2 are shifted 
to higher pt than the maxima in i?AA(A</>). This sug- 
gests that the two effects may not directly related, but we 
observe that the maxima in the Raa{pt) distributions in 
Figs. fT5]fTB1 shift to larger j?x for smaller A<p. To better 
illustrate the shift of the maxima in Raa(pt) we show in 
Fig. [H]the i?AA {pt ) values for the different A<fi bins and 
indicate the variation of the peak position obtained using 
polynomial fits to the first four px bins. For the 30-40% 
centrality bin, the maximum in Raa(pt) for A(f> < tt/12 
is shifted by 0.4 GeV/c relative to the 5vr/12 < A<p < tt/2 
bin. This shift in the peak Raa (pt ) with A</> can pro- 
duce a V2(pt) that peaks higher in px than the inclusive 

RaA- 

The observed shift in the peak of Raa (pt ) with A(f> 
illustrates an important property of collective motion of 
the medium. The collective motion does not simply su- 
perimpose azimuthal variation on the particles produced 
at a given px, it provides a Acj) dependent shift and/or 
broadening in the transverse momentum spectrum of the 
produced particles. The resulting distortion will be the 



smallest for particles produced at angles normal to the 
reaction plane and will be largest for particles produced 
in the plane. Any collective shift of soft particles to 
higher px will increase the measured i?AA(A0,pr) for 
small A(f> relative to large A(j) values producing a simul- 
taneous increase in both the A^-integrated -Raa and the 
v 2 . With increasing px, an expected decrease in the soft 
contamination would naturally explain the simultaneous 
reduction in v 2 and A</>-integrated Raa evident in the 10- 
20% and 20-30% bins where the separation between the 
Raa {pt ) curves for different A</> bins decreases while the 
average Raa also decreases. We will return to investigate 
this correlation again below. 

The 40-50% and 50-60% centrality bins in Fig.[T5]show 
little of the peaking near 2 GeV/c, especially in A<fi bins 
not aligned with the reaction plane. Nonetheless the v 2 
values for the more peripheral bins reach the same large 
maximum values, v 2 ~ 0.2, at intermediate px as the v 2 
values for more central bins where the peak in Raa (pt ) 
is more prominent. Thus, while the peaking in Raa {pt ) 
is less prominent in the more peripheral bins, the relative 
difference between the in-plane and out-of-plane ir° yields 
in the 40-50% centrality bin is comparable to that in the 
20-30% centrality bin. However, it is possible that the 
large A(f> dependence in the more peripheral bins and 
the apparent persistence of that variation to high px in 
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FIG. 13: 7T° Raa versus angle of emission with respect to the reaction plane for 40-50% centrality. Colors/data points as in 
Fig. El 



the 40-50% centrality bin more directly reflects the larger 
spatial anisotropy of the collision zone in more peripheral 
collisions. The question of whether the tt° suppression 
measurements presented here can be understood on the 
basis of geometry and jet quenching will be more fully 
explored in the following section. 

We have observed above that the px and centrality 
dependence of Raa indicate a correlation between inclu- 
sive Raa and i>2 such that the out-of-plane ir° yields vary 
only slowly with px or centrality while the in-plane yields 
approach the out-of-plane yields with increasing px or 
increasing -ZV part . Such a correlation between these two 
seemingly unrelated quantities merely indicates that the 
yields or Raa of 7r°'s measured in-plane and out-of-plane 
more directly reflect the underlying physics responsible 
for the azimuthal variation than the A0-integrated yield 
or Raa and the amplitude of the A</> modulation, V2- 
Indeed, we have argued above that at higher px the cen- 
trality dependence of i?AA(A^) may reflect the geometry 
of jet quenching. At intermediate px, the Raa (pt ) re- 
sults suggest contamination of the in-plane yields by soft 
production and a simultaneous decrease in Raa and V2 
with increasing px as the relative contribution of collec- 
tive soft processes to ir° production decreases. To more 
directly demonstrate the correlation that forms the ba- 
sis of these arguments we show in Fig. [2D] a plot of v-i 



versus the inclusive Raa for centralities from to 60%. 
Data are displayed for px > 2 GeV/c. The intermedi- 
ate centrality bins show a correlated increase of V2 and 
Raa consistent with the discussion above and a possi- 
ble saturation of vi for larger -Raa values. In fact, the 
trends for different centrality bins appear to be in general 
agreement. However, their exact relationship and estab- 
lishing or excluding a causal connection requires further 
investigation. 



C. Nuclear modification factor dependence on path 
length 

The centrality of a collision fixes the geometry of the 
overlap region between the nuclei, and fixing the angle 
of emission of the particles further constrains the path 
length through the medium. We can use this feature to 
study the dependence of the nuclear modification factor 
on the path length traversed by the partons. We investi- 
gate the path length dependence by expanding on several 
methods previously described in [25| . We start with three 
estimators of the path length that are purely geometric, 
and one that includes the color density of the medium in 
its calculation: 

1. We start by modeling the overlap region as an el- 
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FIG. 14: 7T° Raa versus angle of emission with respect to the reaction plane for 50-60% centrality. Colors/data points as in 
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lipse defined by 



? 2 

x z y z 
a z 



b 2 



(14) 



where the minor axis b is oriented in the x direc- 
tion and is parallel to the reaction plane. The 
axes a and b are fixed by the intersection in the 
transverse plane of two hard spheres with R = 
6.8 fm. In terms of the spatial eccentricity e = 
(y 2 — a; 2 ) / (y 2 + x 2 ) (often used with Glauber cal- 
culations), we can express the distance from the 
origin to the edge of the ellipse at a given angle: 



L € (A<t>) 



by/T 



VI + ecos2A0' 



(15) 



Since this length starts at the origin, and does not 
take into account color density, the expression pro- 
vides a very simple estimator with which we can 
evaluate the dependence of the -Raa on path length. 
We will refer to the hard sphere result as L ei hs- 

2. Instead of an ellipse strictly defined by the trans- 
verse profile of two hard spheres, we model the colli- 
sion region as an effective ellipse whose dimensions 
are determined by equating the RMS radius and 



eccentricity to the corresponding quantities calcu- 
lated from the transverse distribution of partici- 
pant density based on standard Glauber calcula- 
tions. This length, L e , is evaluated using the same 
expression as Eq. [T51 with b — yfx 2 . Both quan- 
tities are determined using the PHENIX Glauber 
model [EH. 

3. For a more realistic approach, we evaluate the dis- 
tance along the parton's path weighted by the par- 
ticipant density, 

/•oc 

pL xy = / dlp palt (xo + /cos A(j>, y + lsmA(f>), (16) 
Jo 

where (xo,yo) is the hard-scattering position and 
A(j> is the angle of the jet with respect to the x axis. 
The jet production point is sampled from a Monte 
Carlo using a weighted Taa(^, y) distribution and a 
uniform A(f> distribution. The participant density, 
assumed to be proportional to the color density, 
is calculated from the Glauber model. The yO pa rt 
density in Eq. \W\ is modeled using a ID Bjorken 
expansion, 



p c {r) oc 



7-o 2 



(?) 
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FIG. 15: (Color online) Semi-log plots of Raa{pt) for each A<f> bin, in different centrality ranges. The A<j> bins are represented 
as follows: closed (black) circles, 0-15°; closed (red) squares, 15-30°; closed (light green) triangles, 30-45°; closed inverted (blue) 
triangles, 45-60°; open (magenta) circles, 60-75°; and open (dark green) squares, 75-90°. The systematic error in the inclusive 
-Raa is represented by the grey bands. Errors due to the correction factor have been omitted for clarity. 



Thus pL xy roughly represents LPM energy loss [H[ . 
Note that p c approaches the same 1/r dependence 
as the standard p oc tq/t but differs from by a fac- 
tor of 2 at t = to (additionally this form is regular 
at t = 0). 

4. Finally, we modify pL xy by normalizing it by the 
value of the participant density at the center of 
the collision region, p ccnt = p pa it(0,0). As a re- 
sult, pL xy / p C ent is an estimator based on geometry 
alone, but also accounts for the effect of the den- 
sity distribution both on the jet production point 
as well as the path from that point to the edge of 
the medium. 

Figs. |2"D|2"41 show the Raa dependence on L t: hs, L e , 
pL xy , and pL xy / p cent respectively. The results shown in 
this paper cover the pt range pt = 1-10 GeV/c, not only 
extending the measurement presented previously but al- 
lowing a much finer binning in px ■ The statistical errors 
in the Raa measurements are represented by error bars 
(see Section llV Bp . The systematic errors shown in these 
data are on the Raa values only, and are indicated by 
the filled boxes. The major contribution to the system- 



atic error in the L e values is the determination of -/V par t, 
and is at the 10-20% level. 

Both the L tt hs and L e behavior show an interest- 
ing degree of scaling. This result is all the more un- 
expected because of the overly simple geometric pic- 
ture they represent. Despite the simplistic picture they 
present, they both exhibit striking universality: the hard 
sphere i?AA(£e,hs) scales well in the low pt region (as 
high as pt ~ 4 GeV/c) while i?AA(^e) scales well to 
higher pt, at least one bin in pt beyond L e j ls (though 
one might argue qualitatively this trend continues even 
higher when the most peripheral centrality is excluded). 
The more precise pr dependence available in this data 
set reveals a slight deviation from the universality with 
L c that was previously reported [25| . 

By contrast, we expect the pL xy estimator to provide 
a somewhat more intuitive and concrete picture, as it 
represents a more realistic approach to the geometry and 
medium. Since we expect radiative energy loss to play a 
greater role at high pr, this should be the estimator that 
would provide the best scaling. In fact, at the higher 
Pt range, a universality does emerge, though not until 
Pt ~ 6 GeV/c. Below that value, the measured Raa 
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points lie on parallel, but separate, curves. When pL xy 
is normalized to the central density, data again exhibit 
a more universal dependence over a wider pt range than 
what is seen with pL xy alone, as shown in Fig. 1241 

When considered together these results offer a rich pic- 
ture. At low to moderate px simple geometry may play a 
larger role in determining the final level of Raa than con- 
ventionally thought. At higher px the scaling motivated 
by energy loss (pL xy ) describes the data well. We note 
that there are three possible (and not necessarily exclu- 
sive) interpretations: 1) at low to moderate px the com- 
bined effects of the boost due to expansion and fragmen- 
tation are sensitive primarily to the difference in lengths 
traveled by the partons, and only weakly dependent on 
other parameters 2) we need to restrict the analysis of 
the ir° Raa to px > 5 GeV/c to be in the range where 
fragmentation followed by energy loss dominates, or 3) 
the assumption that energy loss does not depend linearly 
on color density (34[ is incorrect and leads a departure 
from scaling with pL xy at low to moderate px- 



V. SUMMARY AND CONCLUSIONS 

We have presented measurements for the azimuthal 
anisotropy of neutral pions in Au + Au collisions at 
^s NJV = 200 GeV. These measurements include the i>2 
and _Raa(A</>) of 7r°s as a function of transverse momen- 
tum and centrality. The v?, has been measured from 
p T w 1-10 GeV/c in eight centrality bins and four com- 
bined centrality bins. The Raa has been measured in 
the same pt range in six centrality bins. In addition, the 
Raa dependence on effective path length through the 
collision region has been presented for five centralities. 

The general trend seen in the V2(pt) data is an initial 
increase in v 2 with increasing px, peaking of the v% in the 
region of 2 < px < 3 GeV/c, followed by a decrease in the 
V2- We have argued that such a trend implies a transition 
from particle production dominated by soft processes to 
a pt region dominated by hard processes. In order to 
quantify the px and centrality evolution of v 2 , we have fit 
its pt dependence to an empirical expression that allows 
for such a transition. While the statistical precision of the 
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FIG. 17: (Color online) n° -RAA(A jr part ) in reaction-plane bins at fixed pr- The three bins are as follows: closed (black) circles 
are -Raa(0 < A<f> < 15°) (in plane), the closed (red) triangles are the -Raa(30 < Atfi < 45°), and open (blue) squares are 
#aa(85 < A<t> < 90°) (out of plane). 



high-j»T data limits the conclusions we can draw from the 
fits, it is clear that the data support the assumption of a 
decreased dominance of soft processes transitioning to an 
increased dominance of hard process with increasing pr- 
The more peripheral bins suggest that a px-independent 
V2 may be reached by ~ 5 GeV/c. 

The differential probes represented by -Raa(A</>) pro- 
vide a more sophisticated handle on the role of geometry 
in the collision region. For example, we see that in mid- 
centrality collisions the suppression of pions out-of-plane 
is approximately the same as the suppression in more 
central collisions. The data suggest that the interplay 
between the two main effects, namely collective flow and 
jet quenching, may take place not only along the expected 
transition in from soft to hard physics, but perhaps 
also azimuthally, with the quenching effects being domi- 
nant along the direction normal to the reaction plane. To 
further shed light on the transition from soft to hard px 



regions, we have fit the maxima of mid-central Raa (pt ) 
in each A(j) range. Between the in and out of plane direc- 
tions, we observe a shift of 0.4 GeV/c in position of the 
peak of the spectrum. 

In order to further clarify the centrality evolution of the 
azimuthal dependence of the 7r° suppression, we have pre- 
sented the Raa as a function of A^ part in fixed px bins, for 
three directions: along the reaction plane, normal to the 
reaction plane, and midway between. For -ZV part > 100, 
the Raa along the normal to the reaction plane is almost 
constant, a trend seen in most px bins. By contrast, the 
Raa nearest to the reaction plane drops by almost a fac- 
tor of two, converging on the out-of-plane value at the 
highest A'part. This important feature may provide the 
most compelling argument for geometry as the source of 
suppression. Since the path length in the direction nor- 
mal to the reaction plane varies slowly with centrality, 
we would expect the Raa to be nearly flat. Conversely, 
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FIG. 18: (Color online) n° RAA(N pa , T t) in reaction-plane bins at fixed pt, with log scale on the y axis. Colors/data as in Fig. 1171 



the in-plane path length will be sensitive to the degree of 
overlap and strongly influence the observed Raa- Thus 
it would lead to a small suppression in peripheral colli- 
sions while eventually converging on the same value as 
seen normal to the plane as the anisotropy vanishes in 
more central collisions. This effect is further borne out 
in the correlation observed between u 2 and the inclusive 
Raa- However, we have also argued that contamina- 
tion from soft production could produce a similar behav- 
ior and we have no independent indication of how far in 
Pt soft contamination might extend. Nonetheless, under 
both interpretations we can conclude that the Raa for 
pions produced along the normal to the reaction plane 
more directly reflects the physics of quenching. Wc also 
conclude that the correlation between Raa and vi makes 
separate treatment of these quantities disadvantageous. 

We have examined Raa as function of the average 
path length of the parent parton through the overlap re- 
gion in the collision, through the estimators £ e! /j S , L e , 
pL xy / Pcent! and pL xy . Each of the first three quantities 



represents a progressively more sophisticated estimator 
for the distance traveled by the parton, with pL xy at the 
end providing a proxy for LPM energy loss. Compari- 
son of the scaling with these three measures of lengths 
seems to suggest that the pion suppression at low- to 
moderate-pT is mostly dependent upon simply the ge- 
ometric length. The estimator that should in principle 
be the most realistic one, pL xy , exhibits good univer- 
sality at the highest values, suggesting that energy 
loss comparisons should be restricted to the pt range 
Pt > 5 GeV/c. The importance of simple geometry at 
low- to moderate-px is further supported when pL xy is 
normalized by the participant density at the center. This 
normalization effectively makes pL xy a length. These ge- 
ometric descriptions offer a description of the suppression 
both at low and high px regions, clearly showing a transi- 
tion between the ranges. The features seen in the Raa as 
function of path length tie in consistently with the obser- 
vations of a transition in the behavior of the measured 
v%. These two observables, vi and Raa(^), analyzed 
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together provide a valuable set of probes for understand- 
ing the processes governing the suppression of yields in 
Au + Au collisions. 
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FIG. 20: (Color online) ty vi vs. inclusive _Raa. The points 
denote bins in pt as follows: triangles 2 < pt < 2.5 GeV/c, 
inverted triangles 2.5 < pt < 3 GeV/c, circles, 3 < pt < 
3.5 GeV/c; squares, 3.5 < pt < 4 GeV/c; open triangles, 
4.0 < p T < 4.5 GeV/c; diamonds, 4.5 < p T < 5 GeV/c; 
crosses, 5 < pt < 6 GeV/c. Centrality bins are indicated 
by the colors: light blue, 0-10%; black, 10-20%; red, 20-30%; 
green, 30-40%; blue, 40-50%; magenta 50-60%. 
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FIG. 21: (Color online) tt 7?aa versus L e based on the hard sphere calculation. Each panel corresponds to a pt bin. Each 
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FIG. 22: (Color online) 7r° Raa versus L e based on the RMS radius. Colors/data points as in Fig. I21I 
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FIG. 23: (Color online) tt Raa versus pL xy . The units of pL xy are participant x fm. Colors/data points as in Fig. I21I 
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APPENDIX A: DATA TABLES 



TABLE I: ir° v 2 for 0-5% and 5-10% centrality. All errors 
are absolute. 



Centrality 
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TABLE II: w° V2 for other centralities. All errors are absolute. 
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+0 
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0.149 


+0 
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0.012 



0-10% 



50-60% 



10-20% 



0-20% 



20-30% 



30-40% 



40-50% 



20-40% 



40-60% 



0-92% 
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TABLE III: Raa vs. path length for (upper) 1.0 < p T < 1.5 and (lower) 1.5 < p T < 2.0 GeV/c. 



Centrality 


A0 


L t 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.90 


0.514 


0.009 


0.005 




15-30 


3.29 


1.93 


0.481 


0.009 


0.003 




30-45 


3.41 


2.00 


0.443 


0.008 


0.001 


10-20% 


45-60 


3.55 


2.11 


0.400 


0.007 


0.001 




60-75 


3.69 


2.18 


0.358 


0.007 


0.004 




75-90 


3.78 


2.24 


0.337 


0.006 


0.006 




0-15 


2.78 


1.42 


0.565 


0.009 


0.004 




15-30 


2.85 


1.45 


0.533 


0.008 


0.003 




30-45 


2.99 


1.52 


0.477 


0.008 


0.001 


20-30% 


45-60 


3.19 


1.64 


0.421 


0.007 


0.001 




60-75 


3.39 


1.74 


0.371 


0.006 


0.003 




75-90 


3.52 


1.82 


0.346 


0.006 


0.005 




0-15 


2.43 


1.03 


0.634 


0.009 


0.004 




15-30 


2.50 


1.06 


0.599 


0.009 


0.003 




30-45 


2.66 


1.12 


0.526 


0.008 


0.001 


30-40% 


45-60 


2.87 


1.24 


0.459 


0.007 


0.001 




60-75 


3.11 


1.33 


0.397 


0.006 


0.004 




75-90 


3.27 


1.42 


0.360 


0.005 


0.005 




0-15 


2.14 


0.70 


0.721 


0.009 


0.005 




15-30 


2.21 


0.73 


0.669 


0.009 


0.004 




30-45 


2.37 


0.78 


0.596 


0.008 


0.002 


40-50% 


45-60 


2.59 


0.87 


0.505 


0.007 


0.002 




60-75 


2.84 


0.95 


0.424 


0.006 


0.005 




75-90 


3.03 


1.03 


0.379 


0.005 


0.007 




0-15 


1.92 


0.44 


0.798 


0.010 


0.009 




15-30 


1.99 


0.46 


0.751 


0.010 


0.007 




30-45 


2.14 


0.49 


0.665 


0.008 


0.002 


50-60% 


45-60 


2.36 


0.56 


0.562 


0.007 


0.003 




60-75 


2.61 


0.61 


0.485 


0.006 


0.007 




75-90 


2.81 


0.68 


0.435 


0.006 


0.010 


Centrality 


Acj> 


U 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.94 


0.570 


0.007 


0.003 




15-30 


3.29 


1.96 


0.539 


0.007 


0.002 




30-45 


3.41 


2.03 


0.486 


0.006 


0.001 


10-20% 


45-60 


3.55 


2.15 


0.424 


0.009 


0.001 




60-75 


3.69 


2.22 


0.370 


0.005 


0.003 




75-90 


3.78 


2.29 


0.341 


0.005 


0.004 




0-15 


2.78 


1.46 


0.677 


0.008 


0.003 




15-30 


2.85 


1.49 


0.626 


0.008 


0.002 




30-45 


2.99 


1.57 


0.552 


0.007 


0.001 


20-30% 


45-60 


3.19 


1.70 


0.475 


0.006 


0.001 




60-75 


3.39 


1.80 


0.397 


0.005 


0.002 




75-90 


3.52 


1.89 


0.358 


0.004 


0.003 




0-15 


2.43 


1.06 


0.758 


0.009 


0.003 




15-30 


2.50 


1.09 


0.716 


0.008 


0.002 




30-45 


2.66 


1.16 


0.622 


0.007 


0.001 


30-40% 


45-60 


2.87 


1.28 


0.512 


0.006 


0.001 




60-75 


3.11 


1.38 


0.420 


0.005 


0.003 




75-90 


3.27 


1.48 


0.370 


0.004 


0.004 




0-15 


2.14 


0.72 


0.865 


0.010 


0.004 




15-30 


2.21 


0.75 


0.802 


0.009 


0.003 




30-45 


2.37 


0.80 


0.691 


0.008 


0.001 


40-50% 


45-60 


2.59 


0.90 


0.561 


0.006 


0.001 




60-75 


2.84 


0.99 


0.455 


0.005 


0.004 




7^ QO 
i o-yu 


o.uo 


i ns 


u.oyo 


U.UUJ 


DDR 




0-15 


1.92 


0.45 


0.929 


0.010 


0.008 




15-30 


1.99 


0.47 


0.874 


0.010 


0.006 




30-45 


2.14 


0.51 


0.757 


0.009 


0.002 


50-60% 


45-60 


2.36 


0.58 


0.635 


0.007 


0.002 




60-75 


2.61 


0.64 


0.523 


0.006 


0.007 




75-90 


2.81 


0.71 


0.462 


0.005 


0.009 
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TABLE IV: Raa vs. path length for (upper) 2.0 < p T < 2.5 and (lower) 2.5 < p T < 3.0 GeV/c. 



Centrality 


A0 


L t 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.92 


0.562 


0.007 


0.003 




15-30 


3.29 


1.95 


0.534 


0.007 


0.002 




30-45 


3.41 


2.02 


0.470 


0.006 


0.001 


10-20% 


45-60 


3.55 


2.13 


0.413 


0.006 


0.001 




60-75 


3.69 


2.21 


0.354 


0.005 


0.003 




75-90 


3.78 


2.28 


0.329 


0.004 


0.004 




0-15 


2.78 


1.45 


0.681 


0.008 


0.003 




15-30 


2.85 


1.48 


0.633 


0.008 


0.002 




30-45 


2.99 


1.56 


0.547 


0.007 


0.001 


20-30% 


45-60 


3.19 


1.69 


0.452 


0.006 


0.001 




60-75 


3.39 


1.79 


0.375 


0.005 


0.003 




75-90 


3.52 


1.88 


0.325 


0.004 


0.004 




0-15 


2.43 


1.05 


0.752 


0.009 


0.003 




15-30 


2.50 


1.08 


0.691 


0.008 


0.002 




30-45 


2.66 


1.15 


0.597 


0.007 


0.001 


30-40% 


45-60 


2.87 


1.27 


0.490 


0.006 


0.001 




60-75 


3.11 


1.37 


0.392 


0.005 


0.003 




75-90 


3.27 


1.46 


0.339 


0.004 


0.004 




0-15 


2.14 


0.72 


0.869 


0.010 


0.005 




15-30 


2.21 


0.75 


0.813 


0.009 


0.004 




30-45 


2.37 


0.80 


0.691 


0.008 


0.001 


40-50% 


45-60 


2.59 


0.90 


0.550 


0.006 


0.001 




60-75 


2.84 


0.99 


0.443 


0.005 


0.004 




75-90 


3.03 


1.08 


0.375 


0.004 


0.006 




0-15 


1.92 


0.45 


0.946 


0.011 


0.009 




15-30 


1.99 


0.47 


0.870 


0.010 


0.007 




30-45 


2.14 


0.50 


0.757 


0.009 


0.003 


50-60% 


45-60 


2.36 


0.57 


0.606 


0.007 


0.003 




60-75 


2.61 


0.63 


0.482 


0.006 


0.007 




75-90 


2.81 


0.70 


0.413 


0.005 


0.010 


Centrality 


Acj> 


U 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.87 


0.504 


0.007 


0.003 




15-30 


3.29 


1.90 


0.467 


0.007 


0.002 




30-45 


3.41 


1.96 


0.421 


0.006 


0.001 


10-20% 


45-60 


3.55 


2.07 


0.360 


0.005 


0.001 




60-75 


3.69 


2.14 


0.315 


0.005 


0.003 




75-90 


3.78 


2.20 


0.288 


0.004 


0.004 




0-15 


2.78 


1.42 


0.626 


0.009 


0.003 




15-30 


2.85 


1.45 


0.573 


0.008 


0.002 




30-45 


2.99 


1.52 


0.499 


0.007 


0.001 


20-30% 


45-60 


3.19 


1.64 


0.409 


0.006 


0.001 




60-75 


3.39 


1.74 


0.326 


0.005 


0.003 




75-90 


3.52 


1.82 


0.287 


0.004 


0.004 




0-15 


2.43 


1.03 


0.713 


0.009 


0.004 




15-30 


2.50 


1.07 


0.654 


0.009 


0.003 




30-45 


2.66 


1.13 


0.550 


0.007 


0.001 


30-40% 


45-60 


2.87 


1.25 


0.459 


0.006 


0.001 




60-75 


3.11 


1.34 


0.364 


0.005 


0.003 




75-90 


3.27 


1.43 


0.318 


0.004 


0.005 




0-15 


2.14 


0.72 


0.832 


0.011 


0.006 




15-30 


2.21 


0.74 


0.758 


0.010 


0.004 




30-45 


2.37 


0.79 


0.665 


0.009 


0.002 


40-50% 


45-60 


2.59 


0.89 


0.534 


0.007 


0.002 




60-75 


2.84 


0.98 


0.439 


0.006 


0.005 




7^ QO 


o.uo 


i ri7 


u.o / u 


U.UUJ 


u.uu t 




0-15 


1.92 


0.45 


0.916 


0.012 


0.011 




15-30 


1.99 


0.47 


0.861 


0.011 


0.008 




30-45 


2.14 


0.50 


0.738 


0.010 


0.003 


50-60% 


45-60 


2.36 


0.57 


0.594 


0.008 


0.003 




60-75 


2.61 


0.63 


0.480 


0.006 


0.010 




75-90 


2.81 


0.70 


0.409 


0.006 


0.013 
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TABLE V: Raa vs. path length for (upper) 3.0 < p T < 3.5 and (lower) 3.5 < p T < 4.0 GeV/c. 



Centrality 


A0 


L t 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.80 


0.430 


0.008 


0.003 




15-30 


3.29 


1.82 


0.405 


0.007 


0.003 




30-45 


3.41 


1.88 


0.355 


0.006 


0.001 


10-20% 


45-60 


3.55 


1.97 


0.307 


0.006 


0.001 




60-75 


3.69 


2.04 


0.266 


0.005 


0.003 




75-90 


3.78 


2.09 


0.243 


0.004 


0.004 




0-15 


2.78 


1.38 


0.534 


0.009 


0.003 




15-30 


2.85 


1.41 


0.497 


0.008 


0.003 




30-45 


2.99 


1.48 


0.431 


0.007 


0.001 


20-30% 


45-60 


3.19 


1.59 


0.360 


0.006 


0.001 




60-75 


3.39 


1.67 


0.305 


0.005 


0.003 




75-90 


3.52 


1.75 


0.267 


0.005 


0.004 




0-15 


2.43 


1.01 


0.622 


0.010 


0.004 




15-30 


2.50 


1.04 


0.582 


0.009 


0.003 




30-45 


2.66 


1.10 


0.498 


0.008 


0.001 


30-40% 


45-60 


2.87 


1.21 


0.417 


0.007 


0.001 




60-75 


3.11 


1.30 


0.328 


0.005 


0.004 




75-90 


3.27 


1.38 


0.285 


0.005 


0.006 




0-15 


2.14 


0.71 


0.788 


0.012 


0.007 




15-30 


2.21 


0.74 


0.745 


0.012 


0.006 




30-45 


2.37 


0.79 


0.625 


0.010 


0.002 


40-50% 


45-60 


2.59 


0.88 


0.514 


0.008 


0.002 




60-75 


2.84 


0.96 


0.407 


0.007 


0.007 




75-90 


3.03 


1.05 


0.355 


0.006 


0.009 




0-15 


1.92 


0.45 


0.880 


0.015 


0.015 




15-30 


1.99 


0.46 


0.814 


0.014 


0.011 




30-45 


2.14 


0.50 


0.720 


0.012 


0.004 


50-60% 


45-60 


2.36 


0.57 


0.589 


0.010 


0.004 




60-75 


2.61 


0.63 


0.490 


0.008 


0.013 




75-90 


2.81 


0.69 


0.421 


0.007 


0.018 


Centrality 


Acj> 


U 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.76 


0.390 


0.009 


0.004 




15-30 


3.29 


1.78 


0.364 


0.008 


0.003 




30-45 


3.41 


1.84 


0.329 


0.007 


0.001 


10-20% 


45-60 


3.55 


1.93 


0.291 


0.007 


0.001 




60-75 


3.69 


1.99 


0.249 


0.006 


0.004 




75-90 


3.78 


2.04 


0.234 


0.005 


0.005 




0-15 


2.78 


1.36 


0.511 


0.011 


0.004 




15-30 


2.85 


1.39 


0.468 


0.010 


0.003 




30-45 


2.99 


1.46 


0.418 


0.009 


0.001 


20-30% 


45-60 


3.19 


1.57 


0.348 


0.007 


0.001 




60-75 


3.39 


1.65 


0.281 


0.006 


0.004 




75-90 


3.52 


1.72 


0.256 


0.005 


0.006 




0-15 


2.43 


1.01 


0.593 


0.012 


0.006 




15-30 


2.50 


1.04 


0.552 


0.011 


0.004 




30-45 


2.66 


1.10 


0.487 


0.010 


0.002 


30-40% 


45-60 


2.87 


1.21 


0.418 


0.008 


0.002 




60-75 


3.11 


1.29 


0.341 


0.007 


0.005 




75-90 


3.27 


1.38 


0.311 


0.007 


0.007 




0-15 


2.14 


0.70 


0.738 


0.015 


0.010 




15-30 


2.21 


0.73 


0.676 


0.014 


0.007 




30-45 


2.37 


0.78 


0.588 


0.012 


0.003 


40-50% 


45-60 


2.59 


0.87 


0.495 


0.010 


0.003 




60-75 


2.84 


0.95 


0.397 


0.008 


0.009 




7^ QO 
i o-yu 


^ o^ 


i 


u.oou 


u.uuo 


o ni 9 




0-15 


1.92 


0.44 


0.856 


0.019 


0.021 




15-30 


1.99 


0.46 


0.791 


0.017 


0.015 




30-45 


2.14 


0.50 


0.682 


0.015 


0.006 


50-60% 


45-60 


2.36 


0.56 


0.577 


0.013 


0.006 




60-75 


2.61 


0.62 


0.477 


0.011 


0.018 




75-90 


2.81 


0.68 


0.392 


0.009 


0.023 
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TABLE VI: -Raa vs. path length for (upper) 4.0 < p T < 5.0 and (lower) 5.0 < p T < 6.0 GeV/c. 



Centrality 


A0 


L t 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.71 


0.336 


0.010 


0.004 




15-30 


3.29 


1.73 


0.319 


0.009 


0.003 




30-45 


3.41 


1.78 


0.292 


0.009 


0.001 


10-20% 


45-60 


3.55 


1.86 


0.261 


0.008 


0.001 




60-75 


3.69 


1.91 


0.234 


0.007 


0.004 




75-90 


3.78 


1.95 


0.213 


0.006 


0.005 




0-15 


2.78 


1.34 


0.461 


0.012 


0.004 




15-30 


2.85 


1.37 


0.430 


0.012 


0.003 




30-45 


2.99 


1.43 


0.381 


0.010 


0.001 


20-30% 


45-60 


3.19 


1.54 


0.322 


0.009 


0.001 




60-75 


3.39 


1.61 


0.284 


0.008 


0.004 




75-90 


3.52 


1.68 


0.258 


0.007 


0.006 




0-15 


2.43 


0.99 


0.559 


0.015 


0.006 




15-30 


2.50 


1.02 


0.527 


0.014 


0.005 




30-45 


2.66 


1.08 


0.448 


0.012 


0.002 


30-40% 


45-60 


2.87 


1.18 


0.379 


0.010 


0.002 




60-75 


3.11 


1.26 


0.322 


0.009 


0.006 




75-90 


3.27 


1.34 


0.279 


0.008 


0.008 




0-15 


2.14 


0.70 


0.697 


0.020 


0.012 




15-30 


2.21 


0.72 


0.658 


0.019 


0.009 




30-45 


2.37 


0.77 


0.583 


0.016 


0.003 


40-50% 


45-60 


2.59 


0.87 


0.483 


0.014 


0.004 




60-75 


2.84 


0.94 


0.411 


0.012 


0.010 




75-90 


3.03 


1.02 


0.372 


0.011 


0.014 




0-15 


1.92 


0.45 


0.826 


0.026 


0.026 




15-30 


1.99 


0.47 


0.804 


0.025 


0.020 




30-45 


2.14 


0.50 


0.722 


0.022 


0.008 


50-60% 


45-60 


2.36 


0.57 


0.611 


0.019 


0.008 




60-75 


2.61 


0.63 


0.534 


0.017 


0.022 




75-90 


2.81 


0.70 


0.480 


0.0155 


0.030 


Centrality 


Acj> 


U 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.67 


0.301 


0.014 


0.007 




15-30 


3.29 


1.68 


0.294 


0.013 


0.005 




30-45 


3.41 


1.73 


0.267 


0.012 


0.002 


10-20% 


45-60 


3.55 


1.81 


0.246 


0.011 


0.002 




60-75 


3.69 


1.85 


0.218 


0.010 


0.006 




75-90 


3.78 


1.90 


0.199 


0.009 


0.008 




0-15 


2.78 


1.31 


0.420 


0.018 


0.008 




15-30 


2.85 


1.34 


0.382 


0.017 


0.006 




30-45 


2.99 


1.40 


0.347 


0.015 


0.002 


20-30% 


45-60 


3.19 


1.50 


0.301 


0.013 


0.002 




60-75 


3.39 


1.56 


0.260 


0.012 


0.007 




75-90 


3.52 


1.63 


0.240 


0.011 


0.010 




0-15 


2.43 


0.98 


0.492 


0.022 


0.011 




15-30 


2.50 


1.01 


0.469 


0.020 


0.008 




30-45 


2.66 


1.06 


0.428 


0.019 


0.003 


30-40% 


45-60 


2.87 


1.17 


0.370 


0.017 


0.003 




60-75 


3.11 


1.24 


0.315 


0.015 


0.009 




75-90 


3.27 


1.32 


0.318 


0.014 


0.014 




0-15 


2.14 


0.70 


0.750 


0.035 


0.021 




15-30 


2.21 


0.72 


0.708 


0.032 


0.016 




30-45 


2.37 


0.77 


0.557 


0.028 


0.006 


40-50% 


45-60 


2.59 


0.87 


0.466 


0.023 


0.006 




60-75 


2.84 


0.94 


0.385 


0.019 


0.019 




7^ QO 


^ o^ 


1 09 




o oi 7 


097 




0-15 


1.92 


0.45 


0.757 


0.039 


0.050 




15-30 


1.99 


0.47 


0.715 


0.040 


0.036 




30-45 


2.14 


0.50 


0.672 


0.038 


0.013 


50-60% 


45-60 


2.36 


0.57 


0.642 


0.033 


0.014 




60-75 


2.61 


0.63 


0.605 


0.032 


0.038 




75-90 


2.81 


0.69 


0.563 


0.030 


0.051 
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TABLE VII: Raa vs. path length for (upper) 6.0 < p T < 7.0 and (lower) 7.0 < p T < 8.0 GeV/c. 



Centrality 


A0 


L t 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.68 


0.281 


0.017 


0.011 




15-30 


3.29 


1.70 


0.278 


0.017 


0.008 




30-45 


3.41 


1.75 


0.261 


0.017 


0.003 


10-20% 


45-60 


3.55 


1.82 


0.264 


0.018 


0.003 




60-75 


3.69 


1.87 


0.241 


0.015 


0.009 




75-90 


3.78 


1.91 


0.240 


0.015 


0.012 




0-15 


2.78 


1.32 


0.396 


0.026 


0.013 




15-30 


2.85 


1.34 


0.368 


0.023 


0.010 




30-45 


2.99 


1.40 


0.357 


0.023 


0.004 


20-30% 


45-60 


3.19 


1.50 


0.312 


0.021 


0.004 




60-75 


3.39 


1.57 


0.277 


0.018 


0.011 




75-90 


3.52 


1.64 


0.274 


0.017 


0.016 




0-15 


2.43 


0.98 


0.537 


0.035 


0.018 




15-30 


2.50 


1.00 


0.496 


0.030 


0.013 




30-45 


2.66 


1.06 


0.437 


0.029 


0.005 


30-40% 


45-60 


2.87 


1.16 


0.376 


0.026 


0.006 




60-75 


3.11 


1.23 


0.261 


0.018 


0.015 




75-90 


3.27 


1.31 


0.249 


0.019 


0.025 




0-15 


2.14 


0.70 


0.734 


0.050 


0.038 




15-30 


2.21 


0.72 


0.625 


0.043 


0.026 




30-45 


2.37 


0.77 


0.546 


0.044 


0.010 


40-50% 


45-60 


2.59 


0.86 


0.465 


0.036 


0.011 




60-75 


2.84 


0.93 


0.377 


0.029 


0.030 




75-90 


3.03 


1.01 


0.363 


0.027 


0.046 




0-15 


1.92 


0.44 


0.564 


0.051 


0.093 




15-30 


1.99 


0.46 


0.597 


0.059 


0.070 




30-45 


2.14 


0.49 


0.572 


0.054 


0.023 


50-60% 


45-60 


2.36 


0.56 


0.631 


0.049 


0.024 




60-75 


2.61 


0.62 


0.717 


0.063 


0.072 




75-90 


2.81 


0.68 


0.635 


0.056 


0.085 


Centrality 


Acj> 


U 


pL X y 


-Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.67 


0.297 


0.032 


0.026 




15-30 


3.29 


1.69 


0.282 


0.037 


0.019 




30-45 


3.41 


1.74 


0.272 


0.034 


0.007 


10-20% 


45-60 


3.55 


1.82 


0.248 


0.030 


0.007 




60-75 


3.69 


1.86 


0.222 


0.031 


0.020 




75-90 


3.78 


1.91 


0.227 


0.030 


0.030 




0-15 


2.78 


1.29 


0.359 


0.037 


0.025 




15-30 


2.85 


1.31 


0.327 


0.040 


0.018 




30-45 


2.99 


1.37 


0.306 


0.035 


0.007 


20-30% 


45-60 


3.19 


1.47 


0.308 


0.033 


0.007 




60-75 


3.39 


1.53 


0.270 


0.033 


0.020 




75-90 


3.52 


1.59 


0.261 


0.029 


0.027 




0-15 


2.43 


0.99 


0.458 


0.055 


0.042 




15-30 


2.50 


1.02 


0.482 


0.052 


0.033 




30-45 


2.66 


1.08 


0.429 


0.063 


0.012 


30-40% 


45-60 


2.87 


1.18 


0.387 


0.044 


0.012 




60-75 


3.11 


1.26 


0.377 


0.057 


0.034 




75-90 


3.27 


1.34 


0.369 


0.045 


0.048 




0-15 


2.14 


0.69 


0.639 


0.086 


0.091 




15-30 


2.21 


0.71 


0.527 


0.079 


0.058 




30-45 


2.37 


0.75 


0.492 


0.062 


0.022 


40-50% 


45-60 


2.59 


0.84 


0.424 


0.066 


0.022 




60-75 


2.84 


0.91 


0.420 


0.062 


0.069 




7^ QO 
i o-yu 


o.uo 


QQ 




u.uuo 


DQQ 




0-15 


1.92 


0.45 


0.922 


0.161 


0.149 




15-30 


1.99 


0.46 


0.895 


0.150 


0.118 




30-45 


2.14 


0.50 


0.651 


0.106 


0.039 


50-60% 


45-60 


2.36 


0.56 


0.598 


0.078 


0.048 




60-75 


2.61 


0.62 


0.424 


0.089 


0.128 




75-90 


2.81 


0.69 


0.351 


0.042 


0.182 
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TABLE VIII: Raa vs. path length for (upper) 8.0 < p T < 9.0 and (lower) 9.0 < p T < 10.0 GeV/c. 



Centrality 


A0 


L t 


pL X y 


Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.67 


0.296 


0.030 


0.017 




15-30 


3.29 


1.68 


0.280 


0.029 


0.013 




30-45 


3.41 


1.73 


0.255 


0.026 


0.005 


10-20% 


45-60 


3.55 


1.81 


0.258 


0.026 


0.005 




60-75 


3.69 


1.86 


0.220 


0.023 


0.013 




75-90 


3.78 


1.90 


0.219 


0.022 


0.019 




0-15 


2.78 


1.33 


0.415 


0.043 


0.024 




15-30 


2.85 


1.36 


0.356 


0.038 


0.016 




30-45 


2.99 


1.42 


0.359 


0.038 


0.006 


20-30% 


45-60 


3.19 


1.52 


0.329 


0.035 


0.006 




60-75 


3.39 


1.59 


0.288 


0.032 


0.017 




75-90 


3.52 


1.66 


0.320 


0.033 


0.027 




0-15 


2.43 


0.98 


0.387 


0.049 


0.026 




15-30 


2.50 


1.01 


0.543 


0.063 


0.029 




30-45 


2.66 


1.07 


0.476 


0.056 


0.010 


30-40% 


45-60 


2.87 


1.17 


0.396 


0.048 


0.010 




60-75 


3.11 


1.25 


0.318 


0.041 


0.026 




75-90 


3.27 


1.32 


0.308 


0.039 


0.037 




0-15 


2.14 


0.70 


0.724 


0.103 


0.072 




15-30 


2.21 


0.73 


0.658 


0.094 


0.053 




30-45 


2.37 


0.77 


0.614 


0.087 


0.022 


40-50% 


45-60 


2.59 


0.87 


0.535 


0.076 


0.025 




60-75 


2.84 


0.94 


0.382 


0.058 


0.065 




75-90 


3.03 


1.03 


0.310 


0.048 


0.086 




0-15 


1.92 


0.45 


0.536 


0.118 


0.212 




15-30 


1.99 


0.46 


0.682 


0.121 


0.192 




30-45 


2.14 


0.50 


0.661 


0.130 


0.065 


50-60% 


45-60 


2.36 


0.56 


0.508 


0.112 


0.047 




60-75 


2.61 


0.62 


0.791 


0.144 


0.193 




75-90 


2.81 


0.69 


0.658 


0.135 


0.215 


Centrality 


A<f> 


U 


pL X y 


Raa 


Stat Error(abs) 


Sys Error(abs) 




0-15 


3.23 


1.69 


0.397 


0.051 


0.028 




15-30 


3.29 


1.71 


0.327 


0.045 


0.018 




30-45 


3.41 


1.76 


0.223 


0.034 


0.005 


10-20% 


45-60 


3.55 


1.84 


0.249 


0.035 


0.007 




60-75 


3.69 


1.89 


0.203 


0.030 


0.021 




75-90 


3.78 


1.94 


0.213 


0.030 


0.034 




0-15 


2.78 


1.33 


0.340 


0.054 


0.038 




15-30 


2.85 


1.36 


0.319 


0.051 


0.025 




30-45 


2.99 


1.42 


0.316 


0.052 


0.009 


20-30% 


45-60 


3.19 


1.53 


0.374 


0.058 


0.010 




60-75 


3.39 


1.60 


0.392 


0.061 


0.027 




75-90 


3.52 


1.67 


0.343 


0.055 


0.032 




0-15 


2.43 


1.02 


0.637 


0.108 


0.059 




15-30 


2.50 


1.05 


0.588 


0.103 


0.044 




30-45 


2.66 


1.11 


0.519 


0.090 


0.018 


30-40% 


45-60 


2.87 


1.22 


0.491 


0.084 


0.022 




60-75 


3.11 


1.31 


0.331 


0.069 


0.054 




75-90 


3.27 


1.40 


0.274 


0.057 


0.075 




0-15 


2.14 


0.70 


0.428 


0.115 


0.158 




15-30 


2.21 


0.72 


0.612 


0.139 


0.155 




30-45 


2.37 


0.77 


0.316 


0.102 


0.026 


40-50% 


45-60 


2.59 


0.86 


0.530 


0.133 


0.039 




60-75 


2.84 


0.94 


0.678 


0.162 


0.125 




7^ QD 
i j-yu 


o.uo 




u.uuy 


n 1 4^ 


n i/i 

U. 14:0 




0-15 


1.92 


0.45 


1.866 


0.486 


0.498 




15-30 


1.99 


0.47 


1.150 


0.310 


0.273 




30-45 


2.14 


0.51 


0.702 


0.282 


0.091 


50-60% 


45-60 


2.36 


0.58 


0.305 


0.192 


0.077 




60-75 


2.61 


0.64 


0.172 


0.073 


0.369 




75-90 


2.81 


0.71 


-0.005 


-0.002 


-0.812 
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